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ABSTRACT 

We present a study of radiative transfer in dusty, clumpy star- forming regions. A 
scries of self-consistent, three-dimensional, continuum radiative transfer models are 
constructed for a grid of models parameterized by central luminosity, filling factor, 
clump radius, and face-averaged optical depth. The temperature distribution within 
the clouds is studied as a function of this parameterization. Among our results, we 
find that: (a) the effective optical depth in clumpy regions is less than in equivalent 
homogeneous regions of the same average optical depth, leading to a deeper pene- 
tration of heating radiation in clumpy clouds, and temperatures higher by over 60 
per cent; (b)penetration of radiation is driven by the fraction of open sky (FOS) - 
which is a measure of the fraction of solid angle along which no clumps exist; (c) FOS 
increases as clump radius increases and as filling factor decreases; (d) for values of 
FOS > 0.6 — 0.8 the sky is sufficiently "open" that the temperature distribution is 
relatively insensitive to FOS; (e) the physical process by which radiation penetrates is 
preferentially through streaming of radiation between clumps as opposed to diffusion 
through clumps; (f) filling factor always dominates the determination of the tempera- 
ture distribution for large optical depths, and for small clump radii at smaller optical 
depths; (g) at lower face-averaged optical depths, the temperature distribution is most 
sensitive to filling factors of 1 - 10 per cent, in accordance with many observations; (h) 
direct shadowing by clumps can be important for distances approximately one clump 
radius behind a clump. 
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1 INTRODUCTION 

An understanding of the star formation process requires an 
understanding of the underlying density distribution in star 
forming regions. In a direct sense, knowledge of the density 
distribution can help distinguish between different potential 
dynamical scenarios for energy injection, collapse, fragmen- 
tation, and outflow. In a more indirect sense, the density 
distribution significantly affects our ability to infer source 
properties through its influence on thermal balance, chem- 
istry, local emission, and the processing of radiation (radia- 
tive transfer) between the emitting region and the observer. 

The density distribution is important for the dust as 
well as the gas. In particular, the dust forms the domi- 
nant source of opacity to visible and infrared (IR) radia- 
tion, and is the dominant source of IR continuum radiation. 
Perhaps more importantly, the dust dominates thermal bal- 
ance by direct interaction with the radiation field (van de 
Hulst 1949), and through collisions with the gas (Greenberg 
1971; Goldreich & Kwan 1974). As the problems of ther- 



mal balance and radiative transfer are non-local, non-linear 
feedback problems, comparison of detailed models with ob- 
servations remain the best choice of reliably inferring the 
source properties. 

Source geometry is a significant problem in modeling 
star-forming regions. While it is normal to assume some ge- 
ometric symmetry, such restrictions are generally not realis- 
tic. In particular, a wealth of observations (e.g. Migenes et 
al. 1989; Dickman et al. 1990; Falgarone et al. 1991; Cesaroni 
et al. 1991; Marscher et al. 1993; Zhou et al. 1994; Plume 
et al. 1997; Shepherd et al. 1997) show that and fragmenta- 
tion (see also Goldsmith 1996; Tauber 1996 and references 
therein). This is supported by dynamical models (e.g. Tru- 
elove et al. 1998; Marinho & Lepine 2000; Klapp & Sigalotti 
1998; Klessen 1997) which naturally produce clumpy and 
fragmented structures. 

Previous work on radiative transfer in dusty, clumpy 
environments has been undertaken (e.g. Hegmann & Kegel 
2003; Witt & Gordon 1996; Varosi & Dwek 1999; Boisse 
1990). However, much of the previous work has concentrated 
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on scattering (e.g. with application to reflection nebulae), 
the detailed methods of solution, and/or some of the fun- 
damental results such as ability of radiation to penetrate to 
apparently high optical depths. In this paper we use a 3-D 
monte carlo radiative transfer model to extend the previous 
work to a wider and different range of physical conditions. 
In particular, we construct a large grid of models in an ef- 
fort to better delineate and understand the effects of clumpy 
media on radiative transfer. By controlling and varying the 
parameterization of the source, we attempt to disentangle 
some of the underlying physical causes of these effects. 

In section two we describe the model. We discuss the 
general effects of clumping in section three. In section four, 
we introduce the "fraction of open sky" (FOS), and discuss 
the effects of number density, FOS, optical depth, filling fac- 
tor, and clump size on the dust temperature distribution. We 
discuss the effects of shadowing in section five. Finally, we 
draw conclusions in section six. 



2 MODEL 

2.1 Monte carlo model and invariant parameters 

We have constructed detailed, self-consistent, three- 
dimensional radiative transfer models through dust. The 
model utilizes a monte-carlo approach, combined with an 
approximate lambda iteration to ensure true convergence 
even at high optical depths. The model has been tested 
against existing 1-D (Egan, Leung, & Spagna 1988) and 2- 
D (Spagna, Leung, & Egan 1991) codes, and in modeling a 
3-D source (Doty et al. 2005) with good success. Since the 
present study concentrates upon opaque molecular clouds 
where far-infrared radiaton should dominate both for exter- 
nal heating and emission, we ignore the effects of scattering. 
Test cases in both one- and multiple-dimensions show that 
scattering plays only a very small role on the temperature 
distribution within the majority of the sources. 

Based upon the input parameters discussed below, we 
solve for the dust temperature and radiation field at each 
point in the model cloud. The computational volume is 
taken to be cubical of size 1 pc. Each model utilized an 
81x81x81 cell grid, yielding a typical resolution of ~ 3xl0 1B 
cm. The region is taken to be a two-phase medium consist- 
ing of high density clumps, and a low density inter-clump 
medium. The clump/inter-clump density ratio is taken to be 
riciump/nintor-ciump = 100 from observations (e.g. Bergin et 
al. 1996). The external radiation field is taken from inter- 
stellar radiation field (ISRF) compiled by Mathis, Mezger, 
& Panagia (1983). Finally, we adopt the dust opacities in 
column 5 of Table 1 of Ossenkopf & Henning (1994), which 
have been successful at fitting observations of both high- 
mass (e.g., van der Tak et al. 1999, 2000) and low-mass 
(e.g., Evans et al. 2001) regions of star formation. 

2.2 Model parameters 

We adopt a uniformly distributed interclump medium 
interspersed with higher density clumps having 
"-dump Muter- clump = 100. The clumps are randomly 
distributed within the computational volume. The number 
of clumps, and the densities of the clumps and interclump 



Table 1. Model input parameters 



Number 


(L ) 


/ 


r c 

(PC) 


T 


1 





0.01 


0.025 


10 


2 


3 


0.1 


0.05 


100 


3 


300 


0.3 


0.1 


n/a 



medium depend upon the filling factor (/), the clump 
radius (r c ), and the face-averaged optical depth (f). 

The filling factor specifies the fraction of the volume at 
high density. It is given by / = Vciumps/Kotai. We adopt 
filling factors of / = 0.01, 0.1, and 0.3 in accord with obser- 
vations (e.g. Snell et al. 1984; Bergin et al. 1996; Carr 1987). 
When all other parameters are kept fixed, a higher / cor- 
responds to a larger number of clumps, and a more nearly 
continuous dust distribution. Furthermore, due to optical 
depth constraint (see f below), a larger / also corresponds 
to clumps and interclump medium of lower density. 

We make the simplifying assumption that all non- 
overlapping clumps are spheres of radius r c . We choose 
r c — 0.025 pc, 0.05 pc, and 0.1 pc in accord with obser- 
vations (e.g. Carr 1987; Howe et al. 1993). Again, the con- 
straint on the optical depth implies that larger clump radii 
yield smaller densities within the clumps. 

The dust number densities are normalized by averaging 
the optical depth over one entire (81 x 81) face of the compu- 
tational cube. We call this the face-averaged optical depth, 
and denote it by f . This was done to simulate the optical 
depth / column density that might be inferred by a very 
large beam, although it has the same effect as normalizing 
to the total cloud mass. The face-averaged optical depth is 
taken to be f = 10, and 100, in keeping with observations 
from extinction studies (Lada, Alves, & Lada 1999). 

Finally, the models are specified by the strength of 
the internal radiation field, specified by the luminosity of 
the central source, L». We have constructed a grid having 
L* = OLq, 3Lq, and 300Lq to represent a starless core, a 
low-luminosity central object, and a high- luminosity central 
object, respectively. However, in this paper, we restrict our 
report to the starless cores only (L* = 0). The others will 
be presented in a forthcoming paper. 

The ranges of parameters specified above led to a grid 
of 54 models. Each model is numbered by a four-digit integer 
iLlflrcIf- The integers, and their corresponding values are 
given in Tabled for reference. As an example, model 1221 
has L, = 0, / = 0.01, r c = 0.025 pc, and f = 10. 



2.3 Invariance with clump and photon 
randomization 

One concern with "random" clump distribution and monte- 
carlo simulation is the reproducibility of results with differ- 
ent realizations of the clump distribution (i.e. different ini- 
tial seeds in the clump generation), and different photon ray 
paths. As a test, we have considered 9 different initial seeds 
for the clump distributions and photon paths respectively. 

To quantify the temperature distribution, we calculate 
a spherical average temperature, given by 
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Figure 1. Spherical average temperature as a function of radial 
position for varying number of initial photons (ray pays) per cell 
face for the same random clumpy model. 




r/ Rout 

Figure 2. The temperature as a function of position within 
a clumpy (solid line) and non-clumpy (dotted line) low optical 
depth (f = 10) model, along the three principal axes. 



JV JV 

(T(r)) = J2Tim/J2 ni - W 

i i 

Here N is the number of cells a distance between r and 
r + Ar from the center, where Ar is the size of a single cell, 
and Ti and n% are the temperature and density of the cells 
respectively. In keeping with the viewpoint that the high 
density clumps mostly affect the local radiation field, while 
the the low density medium best probes the radiation field, 
the average temperature is taken only over the low-density 
cells. 

The differences in (T(r)) between different clump dis- 
tribution realizations are always much less than IK, and 
are on average less than 0.1K. Based upon this (and direct 
comparison of later analyses) the specific realization of the 
density distribution does not affect our conclusions. 

Likewise, we have modified the number and distribution 
of incident photons to test for sufficient monte carlo cover- 
age. The effect of number of photons on for (T(r)) is shown 
in Fig.0 To be conservative, we adopt 100 photons per cell 
face, at which point differences are less than 0.05K (1 per 
cent). 

Similarly, varying the random distribution of photon 
paths causes deviations of < 0.1K, confirming that the cov- 
erage of the cells by ray-paths is sufficient to concentrate on 
the consequences of clumping. 



3 CLUMPING: GENERAL 

In this section we briefly discuss the general effects of clump- 
ing on the dust temperature distribution. For clarity, we 
directly compare a clumpy, externally heated (L* = 0), low- 
density (/ = 0.1), opaque (f = 10), model having small 
clumps (r c = 0.05 pc) to one with a uniform density distri- 
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Figure 3. The fractional difference in (T(r)) between a clumpy 
and non-clumpy externally heated (L* = 0), low optical depth 
(f = 10) model. 

bution and the same optical depth and external heat source, 
so we can directly measure the effects of clumping. We have 
done similar comparisons for all other models, with similar 
results. 

The temperature distributions along the principal axes 
for a representative clumpy and equivalent homogeneous 
model are shown in Fig. [5] To quantify the differences be- 
tween the temperature distributions, in Fig. [3] we plot the 
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Figure 4. Distribution of temperature (solid lines, left-hand 
scale) and density (dotted lines, right-hand scale) for cuts along 
the x-, y-, and z-axes of model 1221. Other models yield qualita- 
tively identical results. 
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Figure 5. Temperature as a function of angle-averaged atten- 
uation ((e~ T )) for three distances from source center (.Rout/3, 
Rout/2, and 2R ut/3). The solid symbols correspond to the tem- 
perature in the low density interclump medium, while the open 
symbols correspond to the high density clumps. This result is for 
model 1221. Other models are qualitatively identical. 



fractional difference in (T(r)) between these two represen- 
tative models. From these two figures, it is immediately ob- 
vious that the inclusion of clumps - even for the same f - 
changes the temperature structure significantly. In particu- 
lar, the clumpy model experiences higher temperatures by 
up to 50 per cent toward the center, and up to ~ 65 per cent 
at intermediate radii. As a result, we conclude that clump- 
ing itself affects the temperature distribution, even when the 
average source mass or column density is held constant. 

This result confirms the previous finding of others (e.g. 
Hegmann & Kegel 2003; Witt & Gordon 1996; Varosi & 
Dwek 1999) that the effective optical depth in a clumpy 
medium is less than the homogeneous value. It also suggests 
that it is important to understand the way in which the pa- 
rameterization of the clumpy density distribution can affect 
the dust temperature. We address these individually below. 



4 CLUMPING: EFFECTS OF PARAMETERS 
4.1 Number density 

In Fig.0|we plot the temperature (solid lines, left hand scale) 
and dust number density (dotted lines, right hand scale) for 
cuts along the x-, y-, and z-axes. The clumps are signified 
by the higher density regions. As can be seen, the clumps 
are resolved, and appear to be of different sizes as the axes 
do not penetrate all clumps along a diameter. The "wiggles" 
and ~ 0.1 — 0.2 K deviations in the temperature distribution 
are not simply indicative of the uncertainties in the monte 
carlo calculation (see previous results, and Fig. |2J- Instead, 
the majority are dominated by local differences in radiation 
field due to different amounts of blocking of radiation by the 
surrounding clumps (see Sect. 5 for more discussion). 

As can be seen in Fig. 0] the dust temperature is lower 



within the high density clumps than it is within the lower- 
density interclump medium. Since the total emission by 
grains of absorption efficiency (Q^s^a 2 ) oc i/ 13 goes as T 4+l3 , 
one might expect a factor of 100 increase in density to lead 
to a decrease in temperature by a factor 100 1 ' ' 4+ ' 3 '. Given 
that p ~ 1.8 in the FIR for the adopted dust model, this is 
a decrease by a factor of ~ 2.2. For comparison, the average 
decrease within the clumps is a factor of 1.5 with a great- 
est decrease a factor of 1.6. This is due to the fact that the 
clump centers can be warmed by absorbing the radiation 
emitted from the clump edges, in effect trapping radiation 
within the clumps (similar to the opaque centers of centrally 
heated envelopes, Doty & Leung 1994). 

Of even further interest is the fact that while the tem- 
perature decreases coincide with the locations of the clumps, 
the shape of the temperature profiles do not match the steep- 
ness of the clump/interclump interfaces. This further sug- 
gests that radiative transfer effects play a role. 

4.2 Effective optical depth 

Previous authors (e.g., Hegmann & Kegel 2003; Witt & Gor- 
don 1996; Varosi & Dwek 1999) identified the effective opti- 
cal depth (— ln[I/AT/£*]), where Lat is the attenuated flux 
integrated over all solid angles and L* is the unattenuated 
flux integrated over all solid angles or equivalently the aver- 
age attenuation ((e _T ) = J e~ T d£lj J dQ) as a key measure 
of the ability of radiation to penetrate. This is reasonable, 
as the fewer photons penetrating leads to cooler dust, which 
is confirmed in Fig. where we plot the dust temperature 
at three different positions along the x-axis as a function of 
(e" T ) for model 1221. 

One question is left unanswered - the physical process 
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driving the angle-averaged attenuation must be identified. 
There are two possibilities: (1) the ability of radiation to 
penetrate is dominated by general attenuation by multiple 
clumps along a given line of sight and attenuation by the 
low-density interclump medium; or (2) the radiation streams 
primarily through the "holes" between clumps. We consider 
these two possibilities in the following subsection. 



4.3 Optical depth and fraction of open sky 

4- 3.1 Motivation and definition of (r) and FOS 

In order to answer the question of whether and how 
much geometry matters in determining the local temper- 
ature/radiation field, we consider two limiting cases. These 
are a measure with limited geometry information, the angle- 
averaged optical depth ((t)), and a measure which is dom- 
inated by geometry information, the fraction of open sky 
(FOS). We consider these seperately below. 

If the radiation generally diffuses and is attenuated by 
multiple clumps or the low-density interclump medium, we 
might expect to find a dependence of source properties on 
the optical depth averaged over all directions. As a result, 
we define the angle-averaged optical depth to be 



(r) = / r(fi)dn/ / dn 



(2) 



On the other hand, if the radiation penetrates mainly 
by streaming through holes between clumps, then the source 
properties should depend more significantly on the fraction 
of the sky that is uncovered by the clumps. Consequently, 
we also define the fraction of open sky, FOS, to be 



FOS 



qdQ/ / dQ. 



(3) 



Here q = 1 if there is no clump along the given line of sight 
(i.e. if the sky is "open" in that direction), while q = if 
there is a clump along the given line of sight (i.e. if the sky 
is "closed" in that direction). Consequently, FOS=l implies 
an open sky with no clumps, while FOS=0 corresponds to 
a closed sky where all lines of sight are blocked by clumps. 

As a first test of the effect of FOS on temperature, in 
Fig. HJwe plot the temperature (solid lines, left-hand scale) 
and FOS (dotted lines, right-hand scale) for cuts along the 
x-, y-, and z-axes, similar to Fig. [I] There is a strong cor- 
relation of dust temperature with FOS. Interestingly, the 
shape of the temperature distribution correlates with FOS 
to a much higher degree than it anti-correlates with the dust 
number density. This strongly suggests that FOS plays a key 
role, which we test in the following discussion. 



4-3.2 Effect of (t) on temperature 

To test the role of angle- averaged optical depth, we consider 
the dependence of temperature on (r) for a fixed FOS. In 
this way, we can isolate the effects of (r) from FOS. The 
FOS value here was chosen to maximize the number of data 
points to provide for the best possible statistics. The re- 
sulting dependence of temperature on (r) for three different 
radial positions is shown in Fig. [7] We note that the cho- 
sen FOS range is small to make the result independent of 




[R/ Rout] 

Figure 6. Distribution of temperature (solid lines, left-hand 
scale) and FOS (dotted lines, right-hand scale) for cuts along the 
x-, y-, and z-axes of model 1221. Other models yield qualitatively 
identical results. 
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Figure 7. Temperature as a function of angle- averaged opti- 
cal depth ((t)) for three distances from source center (ilout/3, 
Rout/2, and 2R ou t/3), for a fixed FOS. This result is for model 
1221. Other models are qualitatively identical. 



FOS, though moderate increases in the range lead to similar 
results. 

There appears to be little correlation of temperature 
with (r) at a given position. We quantify this by fitting 
the data distribution at each radial distance with a best-fit 
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Figure 8. Slope of the best fit temperature - (r) relationship 
divided by the uncertainty in the slope as a function of position 
in the source. Notice the small correlation that exists, suggesting 
that (t) is of only minimal importance in determining the tem- 
perature distribution. The dotted lines denote 9({t)) = ±3 (i.e. 
a 3<r deviation from zero). This result is for model 1221. Other 
models are qualitatively identical. 



straight line. To understand the significance of the slope / 
correlation, we define a correlation parameter, 



<9(<t}) = m({r))/cr m((T)) . 



(4) 



Here, m((r}) is the slope of the best-fit line relating the 
temperature and (t), and (T m (/ T \) is the uncertainty in that 
slope. The results are shown in Fig.|5] We have also included 
dashed lines at m = ±3<7 m (( T )) as a guide. 

To best understand the utility of this comparison, con- 
sider a case for which the temperature and (r) are uncorre- 
cted. In this case, the fit between them should yield a zero 
slope, and thus 0((t)) = 0. Likewise, a significant correlation 
should yield 9 > 3 (i.e., a 3cr result). The results in Fig. [HI 
are much more consistent with 0({t}) — throughout much 
of the cloud, and only reach 6({t)) = 3 for r > 0.7R ou t- 
Consequently, we conclude that there may exist only a weak 
correlation between the dust temperature and (t). 



4-3.3 Effect of FOS on temperature 

In a similar manner to (r), we consider the effect of FOS 
on temperature. In this case, we keep (r) fixed. In analogy 
with before, (r) was chosen to maximize the number of data 
points to ensure the best possible statistics. The resulting 
dependence of temperature on FOS for three different radial 
positions is shown in Fig. |5] 

Inspection of Fig. ^suggests that there exists some cor- 
relation of temperature with FOS. In a similar manner to 
(r), we quantify this by fitting the data distribution at each 
radial distance with a best-fit straight line. We again define 
a correlation parameter, 
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Figure 9. Temperature as a function of fraction of open sky 
(FOS) for three distances from source center (-Rout/3, -Rout/2, 
and 2iJ out /3), for a fixed (t). This result is for model 1221. Other 
models are qualitatively identical. 
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Figure 10. Slope of the best fit temperature - FOS relationship 
divided by the uncertainty in the slope as a function of position 
in the source. Notice the large correlation that exists, suggesting 
that FOS has significance in determining the temperature distri- 
bution. This result is for model 1221. Other models are qualita- 
tively identical. 
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Table 2. Temperature dependence on FOS and (t). See text for 
discussion. 



Model 


2- max 


*cut«T» 


Zcut(FOS) 


0(FOS)/0«r» 


1111 


0.00 


0.64 


0.77 


n/a 


1121 


0.00 


0.47 


0.84 


n/a 


1131 


0.00 


0.84 


0.64 


n/a 


1211 


0.94 


0.64 


0.30 


2.1 


1221 


0.55 


0.69 


0.27 


CO 


1231 


0.00 


0.74 


0.99 


n/a 


1311 


0.99 


0.10 


0.30 


1.4 


1321 


0.97 


0.97 


0.22 


1.1 


1331 


0.86 


0.94 


0.92 


n/a 


1112 


0.00 


0.64 


0.94 


n/a 


1122 


0.00 


0.49 


0.94 


n/a 


1132 


0.00 


0.82 


0.84 


n/a 


1212 


0.94 


0.74 


0.64 


1.0 


1222 


0.55 


n/a 


0.99 


n/a 


1232 


0.00 


0.79 


n/a 


n/a 


1312 


0.99 


0.07 


0.59 


0.8 


1322 


0.97 


0.97 


0.30 


1.8 


1332 


0.86 


0.94 


n/a 


n/a 



0(FOS) = m(FOS)/a m{FO s). (5) 

Here, m(FOS) is the slope of the best-fit line relating the 
temperature and FOS, and cr m (Fos) is the uncertainty in 
that slope. The results are show in Fig. 1101 We have also 
included dashed lines at ±3cr m ( FO s) as a guide. 

The results in Fig. HOI are generally not consistent with 
6>(FOS) = 0, and lie at or above 0(FOS) = 3 for a good 
deal of the cloud (r > 0.25i? ou t)- Consequently, we can say 
that there exists a significant correlation between the dust 
temperature and FOS. 

Finally, it is encouraging to note that the correlation 
between temperature and FOS is positive. This is expected, 
since a higher FOS leads to more direct heating by the ex- 
ternal radiation field, and is evidenced by m(FOS) > 0. 

4-3.4 Extension and interpretation of (r), FOS results 

We can extend the discussion of FOS and (r) to the remain- 
der of our family of models. The results are summarized in 
Table |21 In this table, we specify the fractional radial posi- 
tion from cloud center as x = r/R out . 

The results in Tabled require an understanding of the 
range over which FOS or (r) can have a significant effect 
on the temperature distribution. In the table, :r max is the 
point beyond which, FOS does not play a significant role on 
the temperature distribution. We find (see Section 4.5) this 
occurs for FOS > 0.6 - 0.8. 

To understand the origin and location of this region, 
consider a point a distance r from the center of the re- 
gion of radius Rout- The number of clumps exterior to r 
is TV = /(i? ou t/r) 3 (l — x). If these clumps are distributed 
about this point on the outward facing half of an equal 
volume sphere of radius r cff = [4i^ ut (l — x) /(An/3)] 1 ' 3 , 
then the fraction of open sky in the outward direction is 
FOS + = 1 - Nr 2 c /2rl s . We plot the resulting estimate of 
FOS + as a function of position in Fig. 1111 As can be seen, 
all models with a low filling factor (xlxx) have FOS + > 0.8, 



f = 0.01 r 0.025 pc 

f = 0.1 a r c = 0.05 pc 



f = 0.3 □ r c = 0.1 pc 




CO 0.4 - / i 

0.3 - / 
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r / R„ut ( = x ) 

Figure 11. An estimate of the FOS in the outward facing hemi- 
sphere as a function of radial position. The solid lines correspond 
to / = 0.01, the dotted lines to / = 0.1, and the dashed lines 
to / = 0.3. Lines without symbols correspond to r c = 0.025 pc, 
while those with triangles correspond to r c = 0.05 pc, and those 
with squares to r c = 0.1 pc. See text for details. 



as does model (x23x) . For these models, FOS is not expected 
to play a significant role in determining the temperature dis- 
tribution. As a result, the radius up to which FOS should 
play a role (= :r max ) becomes zero for these models in Table 
[5] For other models, < x max < 1. 

To understand the comparative importance of FOS vs. 
(t), we consider the average deviation of 6 from zero for 
x < :r max . This is reported in the last column of Tabled For 
models in which a; max = 0, such an average is not possible, 
and the result is noted as "n/a". For the case where only 
#(FOS) > 3 for x < Xmax, the only meaningful correlation 
of temperature is with FOS, and so the result is noted as oo. 
In all other cases, the last column is set to 0(FOS)/0((r)), 
where the bar signifies an average over all positions for which 
9 > 3 and x < x max . As can be seen by the last column 
of Table |21 the effect of FOS exceeds the effects of (r) in 
determining the temperature. 

Finally, one may consider the size of the region over 
which FOS and (r) are important. To do this, in Tablc|5|wc 
specify x cut , defined to be the position beyond which 6 > 3. 
Generally a; cu t((T)) > a; cu t(FOS). At low optical depths, this 
is true whenever x max 7^ 0. For this set, the mean values for 
x C ut{{r)) ~ 0.7, and a; cu t(FOS) ~ 0.4. This suggests that not 
only is FOS (and thus photon streaming) more important 
than (r), but that it is also important deeper into the cloud 
than (r). This is not a suprise, as openings are a significantly 
more efficient method of carrying energy deep into the cloud. 
At larger optical depths, the effect is not as strong due to 
the fact that the less-dense interclump medium is capable 
of meaningful attenuation by itself. 

The analysis above demonstrates that the angle- 
averaged optical depth does not provide sufficient informa- 
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Figure 12. Ratio of local spherical average temperature to the 
value at r = -Rout, as a function of radial position, for a model 
with a clump size of r c = 0.025 pc, and f = 100. The different 
lines correspond to different filling factors: / = 0.01 (solid), / = 
0.1 (dotted), and / = 0.3 (dashed). The results are averaged over 
9 different realizations of the clump distribution, and the error 
bars are the statistical uncertainties from the averaging. 



Table 3. 


Temperature 


contrast {T(r- ln 


))/(Tout) :f = 


r c (pc) 


/ = 0.01 


f = 0.1 


/ = 0.3 


0.025 


0.69 ±0.01 


0.64 ±0.01 


0.56 ±0.01 


0.05 


0.70 ±0.02 


0.67 ±0.07 


0.61 ±0.07 


0.10 


0.70 ±0.04 


0.71 ±0.07 


0.72 ±0.05 



tion for describing a star-forming region. This is seen in the 
fact that the correlation of temperature with FOS is gen- 
erally stronger than the correlation with (r), as well as the 
fact that the region over which FOS is important is gener- 
ally larger than the region over which (r) is important. As 
a result, this suggests that indeed the ability of radiation to 
penetrate further into a clumpy source than a homogeneous 
source is due to the streaming of radiation through holes. 



4.4 Filling factor and clump radius 

In this subsection we investigate the effects of filling fac- 
tor and clump radius on the dust temperature distribu- 
tion. Since / and r c both geometrically parameterizatize 
the clump distribution, we account for the range of re- 
alizations in the actual clump distribution by averaging 
(T(r)) over nine different realizations. Finally, we note that 
in this subsection we present the temperature contrast (= 
(T(r in)} /Tout), which minimizes the effect of absolute den- 
sity scaling between models. 



Table 4. Temperature contrast (T(r in ))/(T ou t> : f = 100 



r c (pc) 


/ = 0.01 


/ = 0.1 


/ = 0.3 


0.025 


0.80 ±0.01 


0.67 ±0.01 


0.58 ±0.01 


0.05 


0.80 ±0.01 


0.69 ± 0.03 


0.55 ±0.02 


0.10 


0.81 ± 0.02 


0.71 ±0.03 


0.57 ±0.05 



4-4-1 Filling factor 

Figure WI\ shows a representative variation in the spherical 
average temperature distribution with filling factor (/) for 
r c = 0.025pc and f = 100. The temperature distributions 
for other models in the grid are qualitatively similar. The 
ratio of inner to outer temperature for the full grid (includ- 
ing uncertainties due to different realizations) are shown in 
Tables 01 and |3 

We see from these data that the role of filling factor 
varies with clump radius. First, the dust temperature de- 
creases with increasing / as expected, since a larger / leads 
to a smaller FOS, and thus less heating available to the lower 
density medium and lower temperatures. 

Second, the temperature distribution for the smallest 
filling factor remains essentially unchanged as the clump size 
changes. This is due to the high FOS throughout, so that 
each point is well-coupled to the external radiation field, 
independent of the size of the clumps. 

Third, the effect of filling factor varies inversely with 
clump radius. This is expected, as FOS increases with in- 
creasing r c , and decreasing /. To see this, consider the sim- 
plifying case of viewing from the center of a sphere, in which 
spherical clumps of radius r c are randomly distributed a dis- 
tance Rout from the center, and in which no clumps overlap. 
This yields FOS = 1 - /i? ut/4r c . For / = 0.01 - 0.3, this 
yields FOS = 0.63 - 1 for r c = 0.1 pc, 0.25 - 1 for r c = 0.05 
pc, and ~ — 1 for r c = 0.025 pc. In reality FOS is influ- 
enced by two further competing factors - clump overlap in 
projection (increases FOS), and clump distribution in dis- 
tance (decreases FOS). The actual variation of FOS with 
model parameters is given in Fig. 1131 The result roughly 
agrees with the estimate above. The lower trend of the ac- 
tual results imply that clump overlap is not a significant 
effect. 

These results are somewhat different for models with 
high optical depth, f = 100. As confirmed in Table 01 the 
temperature distribution at these high optical depths de- 
pends significantly on /, independent of r c . 

This independence of temperature on clump radius for 
large optical depths can can be easily understood by con- 
sidering the optical depth of individual clumps for vari- 
ous models. For the ISRF of Mathis, Mezger, and Pana- 
gia (1983), roughly 1/2 of the total energy is contained in 
A < 25/zm, 2/3 in A < 115/zm, and 3/4 in A < 360/zm. The 
optical depths of a clump at these wavelengths for f = 10 
are, 1.1 < r c i ump (25/jm) < 67, 0.1 < r c i ump (115^m) < 6, 
and 0.01 < r c i ump (360^im) < 0.8. On the other hand, for 
f — 100 the optical depths are ten times higher, namely, 
11 < r c iu m p(25/zm) < 670, 1 < r c i ump (115^m) < 60, and 
0.1 < T c i ump (360^m) < 8. This implies that the clumps for 
f = 10 are opaque to 1/2 of the heating radiation, but trans- 
parent to the remaining 1/2. Consequently, / (and FOS) can 
control roughly 1/2 of the heating radiation. However, for 
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Figure 13. Fraction of open sky at the central point in the 
model cloud (FOS c ) as a function of filling factor for models 
lxxx (L* = 0). The lines correspond to different clump radii: 
r c = 0.025 pc (solid), r c = 0.05 pc (dotted), and r c = 0.1 pc 
(dashed). Note that the results are the same for all optical depths, 
for a given clump distribution. Results for different realizations 
of clump distributions are similar. 
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Figure 14. Ratio of local spherical average temperature to the 
value at r = Rout as a function of radial position, for a model 
with a filling factor of / = 0.1 pc, and f = 100. The different lines 
correspond to different clump radii: r c = 0.025 (solid), r c = 0.05 
(dotted), and r c = 0.1 (dashed). The results are averaged over 9 
different realizations of the clump distribution, and the error bars 
are the statistical uncertainties from the averaging. 



r = 100, the clumps are opaque to 2/3 - 3/4 of the heating 
radiation, meaning that / (and FOS) can control a greater 
amount of the energy that penetrates to a given depth. As 
a result, it is not suprising to realize that as the clumps be- 
come more opaque, filling factor and the ability of radiation 
to stream between the clumps becomes more important. 

4-4-2 Clump radius 

Figure ITTI shows a representative variation in the spherical 
average temperature distribution with clump radius (r c ) for 
/ = 0.1 and f = 100. Again, the temperature distributions 
in the other models are qualitatively similar, and the tem- 
perature contrast results are shown in Tables El&H 

We see from these data that the role of clump radius 
varies with filling factor. In comparison to the case of filling 
factor, the smallest clumps yield the lowest temperatures, 
and the largest clumps the highest temperatures, consistent 
with the FOS findings in Fig. E3 

The independence of temperature with r c for / = 0.01 is 
due to the large FOS. In particular, there exists a large num- 
ber of open lines of sight to any point, making clump vari- 
ation unimportant. Conversely, as / increases, clump size 
becomes more important. This is due to the fact that high 
filling factor models have a commensurately greater num- 
ber of small clumps. As discussed previously, many small 
clumps are more effective at covering the sky than fewer 
large clumps - see Fig. 1131 Consequently, clump radius is an 
important factor in determining the temperature profile for 
high filling factors. 

As before, the results differ for regions of higher opti- 



cal depth, r = 100. In Table 0] we can see the variation in 
the temperature contrast. As discussed previously, the in- 
creasing opacity of the clumps further amplifies the effect 
of filling factor and FOS on the temperature distribution, 
leaving the effect of clump radius as unimportant. 

4-4-3 Relative regions of importance 

By combining the results of the two previous subsections 
(e.g. Fig. ll3l and Tables|3|&|3J , we can infer the regions of im- 
portance for filling factor and clump radius. We first consider 
the case of f = 10. At these moderate optical depths, we see 
that clump radius only becomes important for / > 0.1. In 
the case of / < 0.1, clump radius essentially plays no role 
and is dominated by the fact that the FOS is always high 
for such small filling factors. 

On the other hand, we see that filling factor only be- 
comes important for r c < 0.05 pc. For r c > 0.05 pc, the 
temperature is dominated by the fact that the FOS is rela- 
tively insensitive to /, due to the large number of holes left 
by scenarios with small numbers of large clumps. 

It is interesting that the inferred values of / from ob- 
servations routinely fall in the range of a few per cent (e.g. 
Hogerheijde, Jansen, van Dishoeck 1995; Snell et al. 1984; 
Mundy et al. 1986; Bergin 1996). In light of the results 
above, this may not be a suprise. For smaller filling factors, 
/ and r c do not play a significant role in the temperature 
distribution. On the other hand, for larger filling factors, 
the cloud is closer to homogeneous, and more sensitive to 
clump radius than /. As a result, the temperature distri- 
bution is most sensitive to changes in / for / ~ 0.1. While 
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Figure 15. The deviation of the temperature distribution from 
an equivalent homogeneous model, as specified by 0(/) 2 /</> ma x 
(see text) as a function of the filling factor, for models with r c = 
0.5 pc, and f = 10. 

many interpretations of filling factor are not driven by con- 
tinuum observations, these results are relevant since one of 
the dominant thermal regulators for the gas is collisions with 
the dust - meaning that sensitivity to Td ust near / ~ 0.1 can 
lead to sensitivity in line processes near / ~ 0.1. 

As a test, we have run additional models / = 0.03, and 
/ = 0.2 for the model numbers 1x21. In order to quantify the 
deviation from the homogenous model, we have calculated a 
parameter, <f> 2 = i £((T(r 4 )) - (T hom ( ri ))) 2 /<r{r)l om . Here 
the sum is over radial positions, (T/ lom (r , i)) is the spher- 
ical average temperature for the homogenous model, and 
o"(r)hom is the uncertainty in the temperature distribution 
in the homogenous model. In Fig. 1151 we plot the ratio of 
the 4>{f) 2 to the maximum value, ^ ax as a function of fill- 
ing factor. Notice that, indeed, the deviation from the ho- 
mogeneous models peaks for intermediate values of /. In 
particular, the greatest sensitivity to / occurs in the range 
/ = 0.01 — 0.1 - a resultin accord with the ranges of filling 
factors commonly reported. 

When the optical depth increases to f = 100, we find 
that filling factor plays the dominant role over a much wider 
range. As a result, the sensitivity to filling factor near / ~ 
0.01 — 0.1 for t = 10 is not an effect here. At high optical 
depth, we also find that the clump radius is not important. 
This is due to the fact that / has a larger effect on FOS 
than r c does, and that as optical depth increases, the clumps 
become more opaque, thereby blocking a higher fraction of 
the impinging radiation. 

4.5 Range of importance of FOS 

From the previous discussion, FOS is necessary to de- 
scribe the local radiation field and heating within a clumpy 
medium. Here we attempt to infer the strength of the effect 
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Figure 16. Deviation of temperature for a given point from the 
spherical average, as a function of the FOS for that point. Each 
tenth point is plotted. The grid of models having t=10 (i.e. mod- 
els lxxl) are all plotted. Notice the increased spread for FOS 
> 0.6 — 0.8, and the correlation for FOS < 0.6 or so. 

of FOS on temperature. While this has been addressed in 
passing, here we collect and further distill that information 
to help draw more general conclusions. 

As seen in Section 4.3.4 and Fig. El the greatest effect of 
#(FOS) is at large radii. We find that this is generally true 
across our grid of models. As expected, these points tend to 
have the highest FOS. Interestingly, however, the effect of 
FOS is on the situations where the FOS is relatively small. 
This can be inferred in three ways. 

First, we note that the largest values of #(FOS) occur for 
the models that have the lowest FOS. This can be indirectly 
seen in that the models with the largest #(FOS)/#((r)) in 
Table |5] tend to have the larger filling factors. As seen in 
Figs. Ill) fc ll3l higher / corresponds to lower values of FOS. 

We can roughly quantify this by saying that FOS is im- 
portant for those regions where FOS ~ 0.5 or so. To begin to 
quantify this we can consider the temperature distributions 
and deviations in Tables[3]&|31 together with the FOS values 
in Fig. 1131 In particular, we note that meaningful differences 
in spherical average temperature occur between all models 
at the smallest clump radius, that essentially no variation 
occurs with filling factor for the largest clump radius, and 
that meaningful variation may occur as one changes filling 
factor at the intermediate clump radius. From Fig. 1131 wc 
see that the largest clump radius has values of FOS > 0.5. 
On the other hand, for the smallest clump radius models, 
FOS varies from ~ 0.86 to ~ 0.25 as one changes /. Finally, 
for the intermediate clump radius models, as one changes 
from / = 0.01 to / = 0.1 FOS varies from ~ 0.95 to ~ 0.5, 
with uncertain corresponding change in the temperature dis- 
tribution. But, when the difference between / = 0.01 and 
/ = 0.3 is considered, there is a meaningful temperature 
change. Taken together, this suggests that variation in FOS 
for FOS > 0.5 may not be as important as variations in FOS 
when FOS < 0.5 or so. 

A more direct comparison suggests that FOS always 
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Figure 17. Deviation of temperature for a given point from the 
spherical average, as a function of the FOS for that point. Each 
tenth point is plotted. The grid of models having f=100 (i.e. 
models lxx2) are all plotted. Notice the increased spread for FOS 
> 0.6 — 0.8, and the correlation for FOS < 0.6 or so. 



plays some role, but that a change in FOS becomes increas- 
ingly important in determining the temperature for FOS 
< 0.6. For FOS > 0.6 — 0.8, the spread in temperatures in- 
creases. To see this in Figs. 1161 & 1171 we have plotted the 
fractional deviation of the temperature of a point from the 
spherical average as a function of the FOS at that point 
for all models in our grid. The points plotted cover all radii 
in the models. It is clear by inspection that FOS directly 
correlates with the temperature deviation. Interestingly, the 
correlation is strongest for FOS less than ~ 0.6. On the 
other hand, for FOS > 0.6 — 0.8, there does not appear to 
be a direct relationship between FOS and temperature de- 
viation. However, for the larger values of FOS the spread in 
temperatures is indeed greater. 

Taken together, this agrees well with the idea of stream- 
ing of photons through holes in a clumpy medium. For mod- 
els with a small average FOS, a position with a few extra 
lines of sight can receive significant additional heating. On 
the other hand, when the FOS for a point is high, it is well- 
heated by the external radiation field. However, direct shad- 
owing of a point by a nearby clump, and attenuation by 
the (assumed) tenuous interclump medium can play roles. 
In fact, for the high optical depth models (f = 100), the 
ranges of temperture deviations are smaller than for low op- 
tical depth models, due to the greater attenuation by the 
interclump medium. 



5 SHADOWING 

One further impact of clumps is to directly shadow points 
behind them, leading to lower temperatures. This effect can 
be seen directly in Fig. 2] In this externally heated model, 
points just interior to clumps are decreased in temperature 
relative to equivalent points not directly behind clumps. 
Two important questions arise: (1) over what length scale 



does this effect hold?, and (2) is the cause a decrease in FOS 
for these points? 

To answer the first question, we have considered the 
temperature distribution as a function of distance behind 
clumps. The temperature at these points is consistently 
lower than the spherical average by 5 — 25 per cent. To iden- 
tify the shadowinglengthscale, we determine the distance, d, 
behind the clump at which the deviation is 1/2 the maxi- 
mum. A fit of d of the form d = 8 x r c yields 8 = 1.2 ± 0.4. 
As a result, we infer shadowing is important for lengthscale 
of ~ 1.2 clump radii behind a clump. 

To answer the second question - the role of FOS - we 
can reconsider Fig. |S| which plots the FOS and temperature 
as functions of position along three axes for a representative 
model. Notice how well FOS correlates with temperature, 
especially in the shadowed regions behind clumps. This is 
highly suggestive that FOS is the determining factor in shad- 
owing. More quantitatively, we can consider the fraction of 
sky closed (FCS) by the clump. Using the previous nota- 
tion, the fraction of closed sky is FCS = r 2 /(4d 2 ) = l/4<5 2 . 
For 5 = 1.2, this yields FCS c i ump = 0.17, and a correspond- 
ing FOS of 0.83. This is near the limit of FOS= 0.8 we 
inferred previously above which the sky is sufficiently open 
for the temperature to be insensitive to FOS. As a result, 
we conclude that FOS (and thus actual shadowing) is the 
important mechanism directly behind a clump. 



6 CONCLUSIONS 

We have constructed a grid of three-dimensional continuum 
radiative transfer models for clumpy star-forming regions, 
in order to better delineate and understand the effects of 
clumping on radiative transfer. Based upon this work, we 
find that: 

1. The inclusion of clumps - even for a constant total 
mass / average optical depth - can significantly change the 
temperature distribution within the cloud. These differences 
in temperature can be in excess of 60 per cent, and are due 
to the lower effective optical depth for clumpy media relative 
to equivalent homogeneous media (Sect. 3). 

2. The centers of clumps are warmer than would be ex- 
pected on energy density grounds due to radiation trapping 
(Sect. 4.1). 

3. The temperature distribution is driven by the ability 
of radiation to penetrate, and is thus strongly correlated 
with the angle-averaged attenuation (e _T ) (Sect. 4.2). 

4. While there exists an anti-correlation of temperature 
with density (Sect. 4.1), the correlation with fraction of open 
sky (FOS) is stronger (Sect. 4.2). 

5. We find only a weak correlation of dust temperature 
with angle- averaged optical depth, (r) (Sect. 4.3.2). On the 
other hand, there exists a significant correlation between 
dust temperature and FOS (Sect. 4.3.3). This is interpreted 
as the dominance of streaming of radiation between clumps 
over diffusion through them in determining the radiation 
field (Sect. 4.3.4). 

6. The dependence of radiation penetration on FOS ver- 
sus (r) is robust. The stronger correlation with FOS versus 
(t) extends throughout the grid of models and for different 
realizations of clump distribution (Sect. 4.3.4). 
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7. While (t) may be an effect near the cloud edges, FOS 
is more important deeper into the cloud (Sect. 4.3.4). 

8. At low face-averaged optical depths, f = 10, filling 
factor is more important for small clump radii than large 
clump radii. At large optical depths f = 100, filling factor 
is the dominant effect for all situations (Sect. 4.4.1). 

9. The effects of clump size are only important for the 
largest filling factors (/ = 0.3) and lower optical depths 
f = 10. It is unimportant for lower filling factors or larger 
optical depths, f = 100 (Sect. 4.4.2). 

10. For lower face-average optical depths, f = 10, filling 
factor is most important in determining the temperature 
distribution for / = 0.01 — 0.1, in accordance with most 
observations. For very opaque clouds with f = 100, filling 
factor is important over a larger range (Sect. 4.4.3). 

11. FOS increases as clump radius increases and as fill- 
ing factor decreases (Sect. 4.4). 

12. The variation of temperature with FOS is more sig- 
nificant in cases of small FOS (high filling factor or small 
clump radii), while (r) is relatively unimportant (Sect. 4.5). 

13. For FOS > 0.6-0.8 the sky is sufficiently open that 
there is little dependence of temperature on FOS (Sect. 4.5). 

14. Clumps can directly shadow the regions behind 
them. On average, this regime extends to distance ~ 1.2 
times the clump radius behind the clump, where the clump 
only subtends a small fraction of the sky (Sect. 5). 
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